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We present a model of magnetotransport of inhomogeneous conductors based on an array of 
coupled four-terminal elements. We show that this model generically yields non-saturating mag- 
netoresistance at large fields. We also discuss how this approach simplifies finite-element analysis 
of bulk inhomogeneous semiconductors in complex geometries. We argue that this is an expla- 
nation of the observed non-saturating magnetoresistance in silver chalcogenides and potentially in 
other disordered conductors. Our method may be used to design the magnetoresistive response of 
a microfabricated array. 
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I. INTRODUCTION 

The problem of determining the effective conductivity 
of a classical inhomogeneous medium^ is an old one, but 
a comprehensive theory of the magnetotransport in such 
systems is still far from complete. Moreover, much of the 
early work on inhomogeneous conductors concentrates on 
the zero-magnetic-field case since this is analogous to cal- 
culating the polarization of a random dielectric 2 . For the 
classical problem to be appropriate, the mean free path 
of the charge carriers must be much less than the typical 
length scale of the disorder so that Ohm's law is obeyed 
locally in space: 



E(r) = p(r)j(r) 



(1) 



where j is the current density, E is the electric field and p 
is the resistivity tensor. In the case of a simple conductor 
that possesses a single charge carrier and isotropic inho- 
mogeneities, the resistivity tensor acquires the following 
form in a magnetic field H: 
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(2) 



It is clear that this tensor is characterised by just two pa- 
rameters: the carrier mobility p (since the dimensionlcss 
variable (3 = pH) and the scalar resistivity p$. While 
more complex, anisotropic resistivity tensors have also 
been considered in the inhomogeneous conductor prob- 
lemi^, this paper shall be solely concerned with Eq. J2J. 

For inhomogeneous conductors, the effective resistiv- 
ity p cS is defined by (E) t , = p cS (j) v , where (■■■) v spec- 
ifies an average over volume. In the absence of disorder, 
the magnetoresistance AR/R = [p off (H)~ p e g(0)]/ p cS (0) 
is trivially zero, for arbitrary orientation of the mag- 
netic field. However, in general, the magnetoresistance 
strongly depends on whether it is transverse or longitu- 
dinal, i.e. whether the magnetic field is perpendicular or 
parallel to the current. Furthermore, any nonzero trans- 
verse magnetoresistance must be an even function of field 
due to the rotational symmetry about the current axis. 
Thus, AR/R oc H 2 in the low field limit, which is usually 



defined to be (3 <C 1, but the crossover from the low-field 
to high-field regime is not always obvious in disordered 
systems, as will be discussed later. 



Some recent experiments 
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doped silver chalcogenides, Ag 2 +,5Se and Ag 2 +5Te, have 
added impetus to the investigation of this problem. Both 
silver chalcogenides exhibit a positive, transverse magne- 
toresistance that is a linear function of magnetic field 
throughout the temperature range 4.5K to 300K, with 
no signs of saturation up to fields of 60Ti*i In particu- 
lar, the linearity continues down deep into the low-field 
regime /3 < 1. Such behaviour is not what is seen in con- 
ventional semiconductors, where the resistance increases 
quadratically with increasing magnetic field at low fields 
and, except in very special circumstances, eventually sat- 
urates at fields typically of order lTiSiil, corresponding 
to /3 ~ 1. Since the silver chalcogenides are non-magnetic 
compounds, the origin of the large magnetoresistance is 
unclear, although a quantum theory based on the par- 
tial population of one Landau magnetic band has been 
proposeo-^*^. However, the large range in temperature 
over which the phenomenon occurs suggests that one 
should examine large magnetoresistances resulting from 
classical effects, namely the case where the semiconduc- 
tor is highly inhomogeneous. 

The theoretical study of classically disordered conduc- 
tors may be divided into two separate classes: media 
consisting of two or more distinct phases, separated by 
sharp boundaries, and systems that possess continuously 
variable fluctuations in the conductivity. In the first 
class, solutions for a non-zero magnetic field have been 
derived for an isotropic medium with a low volume frac- 
tion c < 1 of insulating spherical inclusions^SiSi, and 
they give a positive linear magnetoresistance in the high 
field limit (3 3> 1, but the increase of the magnetoresis- 
tance AR/R with field is small, being proportional to 
c. An effective medium method has been used to extend 
this solution to higher volume fractions 2 ^, but this re- 
sult is approximate and it is still only applicable to high 
fields /3 > 1. The class of systems with continuously 
varying conductivity fluctuations has only been studied 
for weak, short-range disorder^S. Using an advection- 
diffusion analogy, the effective magnetoresistance is dc- 
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termined to be AR/R ~ <y 4 / 3 /3 2 / 3 for ^7>1, where the 
disorder width 7 -C 1. 

Thus, the main limitation of the current literature is 
that it is generally restricted to media that only deviate 
slightly from homogeneity, so the increase in magnetore- 
sistance is small and anomalous behaviour only occurs at 
very high magnetic fields. Whilst there is an exact solu- 
tion for the effective magnetoresistance in two dimensions 
that yields a linear magnetoresistance, it is restricted to 
the special case of a two-component media with equal 
proportions of each phasaSi2Si2&. However, this does lend 
credence to the hypothesis that classical disorder is the 
cause of the anomalous magnetoresistance of the silver 
chalcogenides. 

In order to attack the problem of strong inhomo- 
gencities, we previously introduced a two-dimensional 
random resistor network modei2£. We used it to show 
that classical disorder is a possible cause of the anoma- 
lous magnetoresistance of the silver chalcogenides, and 
we raised the possibility of using the networks to con- 
struct magnetic field sensors that operate on princi- 
ples similar to extraordinary magnetoresistance (EMR) 
devices^a. 

In this paper, we investigate the galvanomagnetic 
properties of the random resistor network model in detail. 
This model allows one to study the magnetoresistance of 
an inhomogeneous semiconductor across the whole mag- 
netic field range for a variety of disorder. By considering 
voltages and current paths within the network, as well 
as network magnetoresistances, we demonstrate that our 
resistor network is also capable of simulating macroscopic 
media with complicated boundaries. We use this result 
to include contact effects between resistors within the 
network. 

The paper is organised as follows: Section ITT1 describes 
the random resistor network model and derives expres- 
sions for the network magnetoresistance and Hall resis- 
tance. In Sec. IIIII we use the insight gained from study- 
ing the characteristics of small networks to ascertain the 
symmetries of the network magnetoresistance and, thus, 
determine the condition for which the magnetoresistance 
is non-saturating. Next, we examine larger networks by 
studying the magnetotransport of uniform square net- 
works and random square networks in Sec. IIVI Finally, 
we address the ramifications of contact resistances be- 
tween resistors and boundary effects within the resistor 
network in Sec. El before concluding in Sec. IVII 



II. RESISTOR NETWORK MODEL 

We tackle the inhomogeneous conductor problem by 
discretizing the medium into a random resistor network 
and analysing it numerically. Standard resistor networks, 
where the network unit is taken to be a two-terminal ho- 
mogeneous resistor, are inadequate for simulating current 
flow in a magnetic field, since it does not allow the cur- 
rent to flow perpendicular to the voltage drop across a 



resistor. Thus, a network of two-terminal resistors will 
not faithfully represent Eq. |J2J, which requires the local 
current to make an angle arctan(/3) with the local electric 
field. 
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FIG. 1: The network resistor unit consists of a homogeneous, 
conducting disk with four equally spaced terminals. Currents 
l entering the disk are taken to be positive, while the voltage 
differences v between the terminals are considered positive 
when measured in the clockwise direction. 



The simplest resistor network model that takes account 
of the Hall component is a two-dimensional square lat- 
tice constructed of four-terminal resistors, with a mag- 
netic field applied perpendicular to the network. This is 
sufficient for simulating a transverse magnetoresistance, 
but not a longitudinal one since this requires a three- 
dimensional network. Networks with multi-terminal el- 
ements have also been used to study percolating media 
in a magnetic field^i but, unlike our model, the elements 
were restricted to being either insulators or conductors 
of a set conductivity. 

In principle, the network resistor unit can be of arbi- 
trary geometry but, for simplicity, we take it to be a ho- 
mogeneous circular disk with four current terminals and 
four voltage differences between the terminals, as shown 
in Fig. These currents and voltages are related via a 
4x4 matrix z: 

Vi — Zij Lj (3) 

The coefficients Zy can be determined by solving the 
Laplace equation for the electric potential of a homoge- 
neous, conducting disk, using the currents as boundary 
conditions (see Appendix^). Note that this formulation 
implicitly assumes a uniform injection of current into the 
terminals at all magnetic fields. In practice, a magnetic 
field generally perturbs the current at a boundary be- 
tween two different conductors, but we shall neglect these 
corrections for now, and revisit them later in our discus- 
sion of boundary effects in Sec. IVl 

If the terminals are taken to be equally spaced and the 
angular width ip of the terminal is held fixed (we will 
take ip = 0.14 radians in this paper), then the resistor 
impedance matrix has the form: 

/ a b c d \ 

z= p_( d a b c 

Trt I c d a b I w 

\ b c da) 
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Here, p is the disk scalar resistivity and t is the disk 
thickness, while the matrix elements are dependent on (p 
and 0: a = -g(tp) + f 0, b = g(<p) + \0,c= 0.35 - \ 
and d = —0.35 — \0. In the limit ip — > 0, the func- 
tion g(ip) — > oo so that the disk resistance diverges as 
expected. Like Eq. (JSJ) , the impedance matrix z of each 
resistor is characterized by two independent parameters: 
the mobility p and the quantity s = p/(nt). Note that 
the cyclical permutation of matrix elements is associated 
with any n-terminal resistor that is invariant under rota- 
tions of 27r/n radians. There is also the added constraint 
that J2i v i = 0, so we have a + b + c + d = 0. 
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FIG. 2: Schematic diagram of an TV x M network. One termi- 
nal is grounded to provide a point of reference for the voltages, 
and we can disregard the current at the grounded terminal by 
imposing current conservation. Voltages and currents can be 
classified into 2N — 1 longitudinal components , if and 



2M Hall components V" , I" 



To construct an TV x M random resistor network, we 
connect the disks together (e.g. using perfectly conduct- 
ing wires) and then vary p and s for each resistor. Note 
that we can include positive charge carriers (holes) by 
allowing p to be positive as well as negative, whereas 
previous studies of inhomogencous media have generally 
focused on charge carriers of the same sign. In the con- 
text of real materials such as the silver chalcogcnides, we 
can view the network as representing an array of silver- 
ion clusters and voids within the semiconductor. 

A typical TV x M network is depicted in Fig. [21 One 
can define a network impedance matrix Z so that the 
input voltages Vi = Z^jlj, where I corresponds to the 
input currents. The impedance matrix is determined by 
grounding one terminal to provide a point of reference for 
the voltages, and then using Kirchoff's laws to eliminate 
the current at the grounded terminal as well as eliminat- 
ing the internal currents and voltages within the network. 
By classifying the voltages and currents into 2TV — 1 lon- 
gitudinal components V/* , if and 2M Hall components 
V/ 1 , If, the (2M + 2TV-1) x (2M + 2TV-1) impedance 



matrix Z can be written as 

Z = 



(5) 



To determine the magnetoresistance of an TV x M net- 
work, we set if = and completely ground the left side 
of the longitudinal voltages in Fig. |2| while setting V L on 



the right side to a constant potential U. 
resistance Rnm(H) is then given by: 



The network 



Rnm(H) = 



U 



U 



J2 N I L 



(0) 



where the sum over input currents is performed along the 
ungrounded (right) edge. Similarly, the Hall voltages are 



V H = Z 



HL 



Z 



LLy 



(7) 



If we keep the ratio N/M constant and take the limit 
where N — > oo, then the resistor network should give us 
the galvanomagnetic properties of a real material. Equa- 
tion © is difficult to solve analytically for large networks 
and, in practice, we just use Kirchoff's laws to numeri- 
cally solve for all the currents and voltages in the net- 
work, since this allows us to study the current flow and 
voltage landscape within a network. However, consider- 
able insight can be gained from examining the symme- 
tries of Z. 



III. SMALL NETWORKS AND NETWORK 
SYMMETRIES 

In order to elucidate the basic properties of the resis- 
tor network model, we begin by studying small networks, 
namely 1 X M and TV x 1 networks. The simplest net- 
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FIG. 3: Examples of small resistor networks corresponding to 
(a) a single resistor, (b) 2x1 network and (c) 3 x 1 network, 
where rji = (s;, p t ). 



work is a single resistor, shown in Fig. EI a), and this 
yields AR/R = as expected, because we have assumed 
that the disk is a conventional conductor with no im- 
plicit magnetoresistance. Moreover, this behaviour holds 
for generic lxM networks which are simply equivalent to 
chains of two-terminal resistors. However, TVx 1 networks 
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exhibit non-trivial behaviour since they allow for a plu- 
rality of current paths within the network when N > 1. 
A 2x 1 network of identical resistors (771 = r\i in Fig.OJb)) 
yields AR/R oc 1 while a 3 x 1 network of identical re- 
sistors gives a non-zero magnetoresistance that saturates 
when (3 3> 1. Figure 0] reveals that among uniform N x 1 
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FIG. 4: (Color online) Magnetoresistance of uniform N x 1 
networks with set /x, s, where the saturating curves correspond 
to N = 3, 7 and 11, in order of increasing size. The remaining 
curves represent N = 2, 6 and 10, in order of size again. Inset: 
for odd N the saturation level AR(oo)/R ~ iV 2 



networks there is a general trend for even-TV networks 
to have a non-saturating magnetoresistance and for odd- 
iV networks to possess a saturating magnetoresistance, 
where the saturation level AR(oo)/R scales as TV 2 . We 
see that the differences in AR/R between odd- TV and 
even- TV networks diminish as TV increases, with TV and 
TV + 1 lying on the same curve for sufficiently small f3. 
However, note that the TV x 1 network does not repre- 
sent a well-defined system in the infinite- TV limit since 
R N1 {0)^0. 

Studies of TV x 1 networks consisting of two types of re- 
sistors have also yielded interesting patterns in the high 
field behaviour of AR/R. For the 2x1 network shown in 
Fig-HHb), the magnetoresistance saturates when si/ii ^ 
where the saturation value of AR/R depends on 
(si//i — S2M2), but setting //1 ^ fi 2 and S1//1 = s 2 ^2 
always gives a non-saturating magnetoresistance. An- 
other example is the 3x1 network in Fig. [3Jc) which 
always has a saturating magnetoresistance except when 
si/ii = 2s2^2- Similar patterns occur in networks of 
larger TV. 

The emergence of these symmetries can be understood 
by considering the impedance matrix Z . One can demon- 
strate that it has the form: 



Z = S + f3A 



(8) 



where S and A are symmetric and antisymmetric matri- 
ces, respectively. For chains of resistors, including the 
case where we only have one resistor, S and A are always 



independent of /?, but in general they are only constant in 
the limits where j3 — ► and j3 — ► 00. Note that Z must al- 
ways be a symmetric matrix at zero magnetic field, since 
an antisymmetric matrix implies dissipationless current 
flow. To see this, consider the power of the network 



P = 



V T I = 



I T V 



I T Z T I = I T ZI 



Thus, we have P = if Z T = —Z. Since the presence 
of a magnetic field induces dissipationless flow, it makes 
physical sense to have an antisymmetric matrix attached 
to (3. 

The relevant quantity that is used to derive the mag- 
netoresistance is the (2TV — 1) x (2TV - 1) matrix Z LL . 
From Eq. JSJ, it also has the form Z LL = S LL + j3A LL . 
We can write Z LL = /3Zp so that Zp -> A LL as j3 -> 00. 
Moreover, A LL is an odd antisymmetric matrix for all TV, 
so Zp will possess at least one eigenvalue that approaches 
zero at large fields. 

Now, the sum of the input currents along the right 
edge can be written as 
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where w n and A ra are the nth eigenvector and eigenvalue 
of Zp, respectively. Since (3 appears in the denominator, 
all terms will vanish in the high field limit except for the 
singular terms associated with the eigenvalues approach- 
ing zero. Their behaviour will determine whether the 
magnetoresistance is saturating or non-saturating. If we 
assume that only one eigenvalue Ao approaches zero in 
the high field limit, then we have 

J J 1 \ 
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(10) 



The behaviour of Ao is governed by S / [3 so we must 
have Ao oc 1//3 when (3 ^> 1. Therefore, Ao cancels (3 in 
the denominator and the magnetoresistance is only non- 
saturating if wo.i — > at large fields. This explains 
why special configurations of resistors in the small net- 
works give a non-saturating magnetoresistance. In the 
case where the magnetoresistance is saturating, Ao dom- 
inates the electrical transport and Zj — > wqj. More gen- 
erally, Y)^ wo.i will determine the exact dependence of 
Rnn(H) on field as (3 — * 00. 

To be more concrete, let us consider the simple case 
where 
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(11) 



Then the zero eigenvalue has the (normalized) eigenvec- 
tor woi = (— l) l /y/2N — 1 and in the high field limit we 
obtain 
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(12) 
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where 8(13, TV) is a finite-field correction factor that van- 
ishes as (3 — > oo. Therefore we have the high- field resis- 
tance 



(2N-1 if TV is odd 

Rnn(H) OC < 2N-1 
[S{I3,N) 



if TV is even 



(13) 



The magnetoresistance is non-saturating for even-TV 
networks and it saturates for odd--/V networks, in agree- 
ment with Fig.m The fact that Rnn(H) scales linearly 
with TV for odd-TV networks is in apparent contradic- 
tion with the inset of Fig. but consistency is recovered 
once one notes that R(0) oc 1/TV so that, at high fields, 
Rnn(H)/R(0) oc TV 2 for TV > 1. 

Since Eq. (fH5jl is independent of M, we expect it to be 
valid for all uniform networks. However, in the case of 
large random networks, the situation is more complicated 
since there is generally a distribution of eigenvalues that 
approaches zero at high fields. This must be properly 
treated using sophisticated tools such as random matrix 
theory. 
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FIG. 5: (Color online) Magnetoresistance of iV X N uniform 
networks, where the saturating dotted curves correspond to 
TV = 5, 9 and 13, in order of increasing size, and the remaining 
dotted curves are TV = 10 and 14 as curve size decreases. The 
solid curve represents the straight line AR/R ~ 0.35/3 that 
the dotted curves collapse onto when 1 < f3 < 2JV. Inset: for 
odd N the saturation level AR(oo)/R ~ 1.4TV 



IV. LARGE SQUARE NETWORKS 

To investigate larger networks, we focus on TV x TV 
networks because their zero-field resistance remains finite 
as TV — > oo. Of course, a finite zero-field resistance is 
obtained for any TV x M network provided we keep the 
ratio N/M constant as TV — > oo, but square networks are 
chosen for numerical convenience. 



A. Uniform resistor networks 

The simplest square network is where all the resistors 
are identical and in this case the zero-field resistance is 
constant as system size is increased. In accordance with 
Eq. I|13[) . uniform square networks retain the 'odd-even' 
trend of TV x 1 networks, where odd-TV networks display 
a saturating magnetoresistance and even- TV networks ex- 
hibit a non-saturating one (see Fig.[SJ|. The key difference 
is that the AR/R curves collapse onto a straight line for 
[3 > 1 as TV — > oo, while there are no changes to the 
low field ((3 < 1) behaviour, where AR/R oc [3 2 . This 
non-saturating, linear behaviour resembles the magne- 
toresistance of the silver chalcogcnides and it makes large 
networks candidates for sensors of high magnetic fields. 
Note that AR/R is independent of s since it just appears 
as a scaling of R(H) in uniform networks. From the inset 
of Fig. [31 we see that the magnetoresistance saturation 
of odd networks scales linearly with TV as expected. 

From the point of view of experimentally constructing 
TV x TV uniform networks, it is worth examining the effect 
of adding a finite resistance r at the connections between 
elements. One can mathematically show that this is, in 
fact, equivalent to reducing the angular width ip of the 
terminals. From our numerical simulations, we find that 




FIG. 6: (Color online) Magnetoresistance of 9 x 9 and 8x8 
uniform networks, where there is a constant resistance r at the 
connections between elements. Note that AR/R of the 8x8 
network corresponds to the larger curve at given r. The unit 
in which r is measured, ro, is taken to be R(0) of a uniform 
square network with r = 0. 



it does not change whether the magnetoresistance is sat- 
urating or otherwise, but in Fig HJ we see that the size 
of AR/R at a given field decreases with increasing r. 
Additionally, it changes the field scale such that the di- 
vergence of odd and even curves, as well as the crossover 
from linear to quadratic behaviour, is shifted to higher 
fields. The reduction in AR/R as r increases is not sur- 
prising, because in the limit where r — ► oo, the Ohmic 
dissipation in the network is dominated by the connec- 
tions between disks and we effectively recover a network 
of two-terminal resistors, which has AR/R = 0. 

Another important characterization of uniform net- 
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FIG. 7: (Color online) Hall coefficient Rj 
tions j across a uniform 16 x 16 network. It is dependent on 
H as well as s, and it has the general form Rf = sfifj(P). 



works is the Hall coefficient R H . In terms of network 
parameters (see Fig.|2Jl, we have 
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where it is a function of the position across the network 
j = 1, 2 . . . M. By symmetry, we expect Rf = R^ +1 _j- 
The Hall coefficient converges rapidly with increasing sys- 
tem size, so we will restrict our consideration to uniform 
16 x 16 networks. In Fig [7] we observe that it has the 
general form Rf = spfj((3) so, if we take p = (ne/i) -1 , 
R H is inversely proportional to carrier density n like con- 
ventional semiconductors. In contrast to a conventional 
semiconductor, it is also dependent on f3 at low fields and 
the strength of this dependency increases as we approach 
the network edges j = 1, 16. This already hints that net- 
work boundaries play an important role in the magnc- 
totransport of uniform networks, as will be discussed in 
Sec. E| 



B. Random resistor networks 

To model real inhomogeneous conductors, it is neces- 
sary to consider random resistor networks. In this case, 
we take the distribution of // within the network to be 
Gaussian, with width A/i. Since s is always positive, 
we take s — rj 1 , where r\ also has a Gaussian distribu- 
tion of widt h Arj. We c an then define the width of s to 
be As = \J (rj 4 ) — (rf) , where (• ■ •) is an average over 
the Gaussian distribution. A numerical analysis of ran- 
dom N x N networks produces positive magnetoresis- 
tances that depend on the particular network configura- 
tion for small N and, consequently, exhibit a large range 
in behaviour whose variation increases with increasing 
H. However, this range in behaviour diminishes with 



increasing N, as illustrated in Figure [SJ The distribu- 
tions of magnetoresistance at large field clearly show a 
decreasing distribution width as N increases and the dis- 
tribution becomes evenly spread about the mean for suffi- 
ciently large N. Therefore, the magnetoresistance of the 
infinite random network should be given by the average 
magnetoresistance of finite networks. 
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FIG. 8: (Color online) Distributions of magnetoresistance 
taken from 1000 samples for each network size, where we 
have set (/i) = 0, AfiH = 50, and (rj) — (so that 
As/ (a) = 1A/2). 



Fig. El displays our key results for simulations per- 
formed on 20 x 20 random networks. We find that the 
average AR/R is linearly dependent on field and that 
it is strongly dependent on /i, but it is still insensitive 
to s like in the uniform case. This linear dependence 
can be argued on the grounds that current in a strongly 
disordered medium at large fields is forced to flow per- 
pendicular to the applied voltage a significant proportion 
of the time and therefore contributes the Hall resistance 
p xy oc H to the effective magnetoresistance. The charac- 
ter of the mobility distribution determines the size of the 
relative magnetoresistance because, at sufficiently large 
magnetic fields, we see that AR/R oc (p) for A/i/ (/i) < 1 
and AR/R oc A/i for A/x/(/x) > 1, where the exact pro- 
portionality constants depend on the details of the dis- 
tribution. Therefore, we would expect AR/R of an in- 
homogeneous semiconductor to diminish with increasing 
temperature, since this corresponds to a decrease in /i 
due to phonon excitations. This is consistent with ex- 
periments on the silver chalcogcnides4. 

The crossover from linear to quadratic behaviour oc- 
curs at field (p,)^ 1 for A/z/(/z) < 1 and (A/x) -1 for 
A/z/(/i) > 1. Thus, even when the characteristic field 
(/i) _1 is of order IT, the measured crossover field of a 
disordered semiconductor can be several orders of mag- 
nitude smaller, provided A/i is large. This yields a possi- 
ble explanation for why the linearity of the silver chalco- 
genide magnetoresistance continues down to fields as low 
as lOOe. 

It is also of interest to determine how disorder affects 
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the boundary conditions, then the magnetoresistance of 
the uniform network becomes equivalent to that of a sin- 
gle resistor, and is thus zero. So why is the magnetoresis- 
tance that we calculated for the infinite uniform network 
non-zero and non-saturating? The answer lies in bound- 
ary effects due to the perfectly conducting electrodes that 
are used to apply the potential difference across the net- 
work. 
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FIG. 9: (Color online) Average magnetoresistance AR/R, as 
a function of dimensionless magnetic field H/Ho, of 20 x 20 
random resistor networks for 3 different mobility distribu- 
tions, where Hq — IkOe is a typical field scale. The mag- 
netoresistance was averaged over 10 random network configu- 
rations and, in order of increasing size, the curves correspond 
to: (i) (fi) = 0.1H- 1 , Ap =_H~\ (ii) (fi) = Hq 1 , Ap = 0, 
and (iii) (p) = 0, Ap = 5H . Inset: By scaling the curves 
so that they all have the same magnetoresistance at around 
4:Hq, it can be seen that linearity continues down to lower 
fields when the mobility disorder is large, Ap Hq 1 



the Hall coefficient, because experiments on the silver 
chalcogcnides have established that an anomalous uni- 
versal relationship exists between the magnetoresistance 
and the Hall resistance^. Unfortunately, an enormous 
range in behaviour is displayed for the Hall resistance of 
finite random networks, and there is no obvious conver- 
gence in the behaviour for network sizes ./V < 30. Thus, 
we need to examine even larger networks in order to de- 
termine the Hall resistance for the infinite network. One 
possible approach is to implement a numerical renormal- 
ization group technique where each resistor unit is re- 
placed by a new, renormalised resistor unit consisting of 
a 2 x 2 resistor network, but this is beyond the scope of 
this paper. 



V. BOUNDARY EFFECTS 

A. Effects of macroscopic boundaries 

Before we conclude our study of large resistor net- 
works, we need to address an apparent conundrum: the 
uniform network in the infinite limit should behave like a 
classical homogeneous conductor with no magnetoresis- 
tance. To see this, consider a single resistor, like Fig. ^ 
within the infinite uniform network. From translational 
symmetry, current entering the resistor from the right 
(bottom) terminal is equal to the current leaving from 
the left (top). If we assume that the current flowing per- 
pendicular to the applied voltage is zero, as dictated by 
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FIG. 10: (Color online) An L\ x Li homogeneous medium 
with two perfectly conducting electrodes attached. The elec- 
trode boundary may be treated as an interface between an 
ideal conductor and a material with finite resistivity pi. In 
a magnetic field, current j enters the ideal conductor at an 
angle 9 = arctan j3 with respect to the electric field E. 



Fig. ^| depicts the boundary between the ideal elec- 
trode and a material of finite resistivity p\ in the x-y 
plane. If a magnetic field H is applied in the z-direction, 
the classical electrical transport in homogeneous mate- 
rial obeys Ohm's law, with a resistivity tensor p given by 
Eq. J2J). Now, the electric field E y that is parallel to the 
surface must be continuous across the boundary accord- 
ing to the standard Maxwell equations. Since the electric 
field inside an ideal conductor is always zero, then E y = 
and the electric field outside the conductor must there- 
fore be perpendicular to the perfectly conducting surface. 
This, combined with the form of p, causes the current to 
enter and exit the perfectly conducting electrodes at the 
angle 9 = arctan/?. For strong fields j3 3> 1, the current 
is angled at almost 90° to the electric field, so the effec- 
tive resistivity of the material close to the electrodes is 



PeS 



3y 



Pi0 



(15) 



This provides an explanation for the linear magnctorc- 
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FIG. 11: (Color online) Visualizations of currents and voltages at large field in a 20 x 20 network of disks with radii 1 (arbitrary 
units), where the potential difference U = —IV. The black arrows represent currents, where arrow size corresponds to the 
magnitude of the current, (a) Uniform network at /3 = 30. (b) Random network with (fi) = 0, AfiH = 30 and As/ (s) = l/y/2. 



sistance of the infinite uniform network in Fig. |SJ In 
general, currents will be perturbed at a boundary that is 
perpendicular to the x-y plane if there is a mismatch of 
Hall electric fields E y across the interface when j y = 0. 

Figure lTTT al demonstrates that, in a large uniform net- 
work, the current is most strongly perturbed at the elec- 
trode boundaries in a strong magnetic field. Deep within 
the network, away from the boundaries, the current is 
uniformly spread, so a measurement of the bulk magne- 
toresistance using a four- probed measurement, will yield 
a zero magnetorcsistancc in the infinite size limit. There 
is also the extra restriction that no current is allowed 
to enter or leave the top and bottom edges of the net- 
work, so this forces the majority of the current to enter 
the network at the top left corner and leave at the bot- 
tom right corner. As (3 — > oo, we obtain singularities of 
the current in the aforementioned corners. This type of 
behaviour has already been noted in the context of real 
homogeneous materials^!. 

In addition, the anomalous Hall coefficient in Fig. 
qualitatively matches calculations for the response of sim- 
ple Hall devices constructed from homogeneous materials 
with a square geometrja^SiS like that in Fig. ^] There, 
a geometrical correction factor is used to describe the 
diminution of the Hall voltage due to finite size effects, 
and this factor is dependent on the ratio of the electrode 
width L\ with respect to the length L2 of the Hall device, 
i.e. the ratio N/M of an infinite uniform resistor network. 

It is important to stress that random resistor networks 
correspond to an entirely different class of system from 
the homogeneous conductor: the infinite random resistor 
network represents an inhomogeneous material, while a 
uniform network may only be regarded as inhomogeneous 
when the system size is finite. Comparing Figs. Illf a) 
and (b), we see that the current paths are highly inho- 



mogeneous and filamentary within the random network 
unlike the uniform case. The voltage landscape is also 
nontrivial and the current paths create loops within the 
random system. Therefore, the magnetorcsistancc should 
be non-zero deep within the random network, away from 
the boundaries. 

We can strengthen this claim by considering the volt- 
age correlation function: 

V corr (x.) = (V(r + x)-V(r)) r (16) 

where V^r) is the voltage at position r in the network, 
x denotes vectors oriented in the x-direction, and (• • • ) r 
represents an average over r. Taking the radii of the 
disks to be 1 (arbitrary units), we should note that 
V corr (Nx) — U in an N x N network, since this cor- 
responds to measuring the potential difference across the 
whole network. 

Figure^lplots the voltage correlation function at large 
magnetic field of 20 x 20 networks for four different cases 
of disorder. When the network is uniform, we see that 
the slope of V corr (x) suddenly increases for large x, which 
implies that the voltage drops within the network are 
smaller than those close to the boundary. This demon- 
strates that much of the magnetoresistance of the net- 
work is confined to the boundaries, as expected. In con- 
trast, as we increase the disorder in the network, the 
behaviour of V r corr (x) tends to a straight line, indicating 
that the linear magnetoresistance is spread across the 
whole of the network and is not just a boundary effect. 
Of particular interest is the fact that maximum insensi- 
tivity to the boundary is achieved when (fx) = 0, which 
could imply that the magnetoresistance is largest when 
electrons and holes are present in equal proportions, as 
has been measured in experiment^. Recent calculations 
by Guttal and Stroud^ on two-dimensional, two-phase 
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FIG. 12: (Color online) Voltage correlation function at large 
magnetic field for 20 x 20 networks of varying disorder. Note 
that the cases with (77} = correspond to As = (s) / v2, since 
s = rj 1 . For random networks, the function has been averaged 
over 100 samples. 



media further support these observations. They prove 
that the magnetoresistance of their inhomogeneous con- 
ductor is linear when (p) = 0, but it can saturate when 
(p) ^ 0, like in experiment. 



B. Finite-element modelling of macroscopic media 

In addition to studies of random media, there is also in- 
terest in understanding classical macroscopic media with 
complex geometries, because of possible geometrical ef- 
fects in magnetoresistance and Hall measurements. Ge- 
ometric enhancements of the magnetoresistance are al- 
ready the basis of sensitive EMR magnetic-field sen- 
sors2&2&. In general, it is difficult to analytically cal- 
culate the magnetotransport of macroscopic media with 
complicated boundaries, because the calculation involves 
solving differential equations for the currents/voltages 
where the boundary conditions contain derivatives that 
arc oblique to the boundary surfaces^*. Thus, standard 
mathematical techniques, such as separation of variables, 
will typically fail in these problems, although the applica- 
tion of conformal mappings in two dimensions has proved 
successful in dealing with simple geometries^&SiSi. How- 
ever, one can, in principle, use infinite uniform net- 
works to simulate two-dimensional, macroscopic, com- 
posite conductors, in a manner analogous to the finite- 
element modelling of EMR devices^!. Moreover, cur- 
rent perturbations at the connections between disks can 
be disregarded entirely in these networks if we choose 
the terminals to have the same resistivity as the disks. 
Boundaries within the macroscopic system still present a 
potential problem since they involve disks of differing re- 
sistivity connected together. But, at fixed magnetic field, 
the magnitude of these contact effects will tend to zero as 
the granularity of the network goes to zero (TV, M — > 00), 



provided the number of elements with contact effects 
scales slower than the total number of elements. This 
is certainly true for one-dimensional boundaries within 
a two-dimensional homogeneous medium, since the num- 
ber of boundary elements in the network scales like ./V 
while the total number of elements scales like TV 2 . 

It is important to note that the homogeneous conduc- 
tor constructed from the infinite uniform network will 
possess a mobility p* and resistivity p* that is different 
from those of the elements that generate the network. 
Generally, these effective network parameters will depend 
on the geometry of the element as well as the connections 
between elements. For the situation where there is no 
resistance between the elements, the effective quantities 
can be determined from the resistivity and Hall coeffi- 
cient of a single element. Therefore, using Eq. we 
have parameters 



P* = ^(.9M + 0.35) 



7T/i 



2(g(tp) + 0.35) 



(17) 
(18) 



The high-field magnetoresistance of the uniform square 
network in Fig. is then given by AR/R ~ p*H. 

To confirm the validity of our numerical approach, 
we can compare the magnetoresistance of infinite uni- 
form networks with results of the conformal mapping ap- 
proach. Following the method of Rendell and GirvinS, 
we find that the L\ x Li homogeneous medium in Fig. llUI 
has resistance: 



R = pi v/l + P 



J dx cos &l 1 l 2 (x) 
Jo dx cos Ql 2 l x (x) 



(19) 



where 
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LiL, 



i(odd) 



4 arctan(/3) sin(n7ra:) 
^ cosh(^) 



Thus, for the special case where L\ = L2, we have the 
exact result R = p\ y 1 + /3 2 , which is in agreement with 
our numerical simulations of TV x TV networks in Fig. [S] if 
we take ./V — > 00, f3 = p*H and p\ = p* . 

For a general TV x M uniform network with M > N, 
the numerical simulations in Fig. 1131 demonstrate that 
the resistance is approximately given by the expression: 



Rnm{H) ~ Rnn(H) 



Thus, in the limit where TV, M — > 00, we have resistance 



M 



R NM ~ p* (y 1 + {p*Hf - 1 + —J (21) 

This becomes independent of network dimensions as 
H — > 00, like previous studies of two-terminal devices 
have predicted^, but the high-field magnetoresistance is 
AR/R oc TV/M since R NM (0) = p*M/N. 
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FIG. 13: (Color online) Comparison of resistances from a rep- 
resentative sample of N x M networks with M > N, where 
we have defined ARmm = Rnm{H) — Rnm(0) and we have 
set s = 1. These curves approach Eq. 12 U in the limit of infi- 
nite network size. Curves with M > 3iV are indistinguishable 
from the M = 3N curves and are therefore not shown. 



C. Contact resistance between elements 

The insights gained in the previous section can be used 
to analyse contact effects between elements in the net- 
work. We previously assumed a uniform current injection 
into the disk terminals, but this is generally difficult to 
achieve when (3 3> 1 . The assumption of uniform current 
is only valid for all fields when the Hall resistance pp and 
thickness t is identical for each clement. To assess the 
ramifications of current perturbations at the disk termi- 
nals, we consider the simple case of ideal metal bridges 
connecting the disks. Our numerical simulations show 
that current distortions are restricted to the vicinity of 
the perfectly conducting electrode, while other calcula- 
tions^ demonstrate that the lengthscale of the current 
distortion is proportional to tp. Therefore, when tp -C 1, it 
is legitimate to replace the bridge with a two-terminal re- 
sistor possessing a field-dependent resistance. Using our 
results for a homogeneous conductor with boundaries, we 
determine the contact resistance between disks i and j to 
be 

= l(^h( IM H) + £h{i* j n)) (22) 

where we take h((3) = \/l + f3 2 — 1 to obtain the correct 
low-field and high-field limits. 

The contact resistance cannot be eliminated from our 
network model by reducing the terminal width (p, because 
the disk impedance coefficients associated with field in 
Eq. are independent of tp like p c , even though the 
disk resistance tends to infinity as tp — > 0. Moreover, in 
the case of the infinite uniform network at large fields, the 
contact resistance p c is equal to the network resistance 
ARnm without metallic bridges for all tp <C 1. 



From simulations of random square networks that in- 
clude this contact resistance, we find that the size of the 
network magnetoresistance is increased by up to 100% 
or more, depending on the network disorder, but the ma- 
jor results are qualitatively unchanged: the variation in 
magnetoresistance decreases with increasing network size 
like in Fig. [SJ the average magnetoresistance is linearly 
dependent on field, and the crossover point is determined 
by the mobility distribution, as in Fig. [5] One impor- 
tant consequence of these contact resistances is that the 
bulk magnetoresistance is always non-saturating and lin- 
ear, plus the network's sensitivity to the network bound- 
aries is reduced. For example, an infinite uniform square 
network with metal bridges will possess a non-zero bulk 
magnetoresistance of p c /R(0) due to the contacts, and we 
find that this accounts for about 70% of a two-terminal 
measurement of the network's magnetoresistance. 

Note that contacts of perfectly-conducting wires repre- 
sent an extreme limit where the contact effects are great- 
est. A potentially richer case is where the wires are re- 
placed by interfaces between elements. Here, the mag- 
nitude of each contact resistance in a random network 
ranges from Eq. (1221) right down to zero, when neighbour- 
ing elements have the same Hall resistance and thickness. 
Thus, we can expect to recover bulk magnetoresistances, 
similar to those displayed in Fig. ^] that sensitively de- 
pend on network disorder. 

The situation is further complicated when we consider 
three-dimensional effects at the connections between re- 
sistors. An interesting example is where the disk resistiv- 
ity p and mobility p are constant within the network, but 
s is varied by altering the disk thickness t. Calculations 
by Bruls et al42iii that involve mapping sample thickness 
variations onto a two-dimensional problem, have shown 
that sharp changes in thicknesses, like those at the inter- 
faces between resistors, will have resistance 

p* ~ P PA tJ (23) 

where Ay = *~ 3 , U > tj and /3Ay ^> 1. Thus, the 
contact resistance still has a linear field dependence at 
large fields, but the size of the effect is reduced so that 
the assumption of uniform current injection is now valid 
for (3 <C A" 1 . A thorough analysis of these more complex 
contact effects will be the subject of future work. 

VI. CONCLUSION 

In this paper we have modelled an inhomogeneous con- 
ductor using a two-dimensional random resistor network 
that consists of /crar-terminal resistors in order to take 
account of the Hall component. We have shown that the 
network impedance matrix Z becomes an odd, antisym- 
metric matrix at large magnetic field, so that the high- 
field behaviour of the magnetoresistance is determined by 
the zero eigenvalue of Z. We find that a non-saturating 
magnetoresistance can be obtained in networks as small 
as 2 x 1, where a plurality of current paths is allowed 
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within the network, while large networks typically pos- 
sess a linear magnetoresistance. This is in contrast to 
EMR devices that exhibit an extremely large but satu- 
rating magnetoresistanc e) 28 ! 29 . 

By considering large square networks, we have demon- 
strated that uniform networks in the limit of infinite size 
are equivalent to homogeneous conductors and the ob- 
served linear magnetoresistance in this system results 
from boundary effects at the macroscopic electrodes. As 
such, they can be used to model macroscopic media with 
complex boundaries. However, large random networks 
model strongly inhomogeneous semiconductors and their 
magnetoresistance is not simply a boundary effect. They 
correctly reproduce the anomalous magnetoresistance of 
the silver chalcogenides: non-saturating behaviour with 
a linearity that continues down to low fields for large mo- 
bility disorder. Moreover, the magnetoresistance may be 
large when the Hall resistance is zero, like in experiment. 

The advantage of such a phcnomcnological model of 
positive, non-saturating magnetoresistance is that it is 
potentially relevant to a whole range of materials. Al- 
ready, similar magnctorcsistances have been observed in 
metallic VCv thin films^, micro-sized Co x -Ci_ x compos- 
ites^ and LaSb2 crystals^. Two-dimensional electron 
gases also show a mysterious linear magnetoresistance^ 
and classical disorder has been cited as a possible caused. 

A major limitation of our random resistor network 
model is that it is restricted to two dimensions and, 
thus, cannot describe longitudinal magnetoresistance. It 
is also known that weakly disordered systems with con- 
tinuous fluctuations in the conductivity possess magnc- 
torcsistances that depend on the dimensionality. There- 
fore, we must extend our resistor network model to three 
dimensions in order to fully simulate an inhomogeneous 
semiconductor. However, we anticipate that our two- 
dimensional resistor networks will motivate experiments 
on the magnetotransport of systems with controlled dis- 
order and on high-field magnetic sensors. 

Finally, it would be interesting to explore the magnc- 
tothermopowcr of our resistor networks. A giant mag- 
netothermopower is associated with Ag2_,5Te samples^ 
and its origin may also lie in classical disorder, because 
experiments on composite semiconductor-metal struc- 
tures demonstrate that the magnctothcrmopower can be 
geometrically enhanced^. 
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APPENDIX A: IMPEDANCE MATRIX OF A 
4-TERMINAL CIRCULAR DISK 

We begin by writing the electric field in terms of the 
potential, E = — VZ7, and then combining charge conser- 
vation V. j = with Ohm's law to obtain the following 
differential equation for the potential: 

d f dU_\ 

dxi \ lk dx k ) 

where the conductivity tensor a = pr 1 . For the case of 
a homogeneous medium, Eq. IjAlf) is simply the Laplace 
equation. 

Consider the homogeneous disk in Fig. ^ If we assume 
uniform current injection into the terminals, then we can 
use the currents as the boundary conditions to solve the 
Laplace equation for the potential. In the absence of a 
magnetic field, it is sufficient to take ip -C 1 in order for 
this assumption to be valid, but the currents will gener- 
ally be distorted when f3 3> 1. To simplify the problem, 
we shall initially neglect these distortions. 

Taking the currents entering each terminal to be Ii, 
I2, I3, and J4, respectively, we then obtain the following 
potential along the edge of the disk: 

00 

U(J3,6) = -^J2—[(S-0T)coB(jtf) 

TTipt * — ' n z 
V 71=1 

+(T + pS) sm{n6)\ (A2) 

where 9 defines the angular position on the disk edge, 
and we have 

S = 2I x $m{mp/2) 

+ /2[sin(n7r/2 + nip/2) — sin(n?r/2 — nip/2)] 

+ /3[sin(n7r + nip/ 2) — sin(n7r — ntp/2)] 

+ /4[sin(3n7r/2 + ntp/2) — sin(3n,7r/2 — nip/ 2)} 

T = I 2 [cos(nir/2 - ntp/2) - cos(mr/2 + n<p/2)] 

+ J4[cos(3ri7r/2 — mp/2) — cos(3n7r/2 + mp/2)] 

To determine the impedance matrix z, we take the po- 
tential differences between the equally-spaced terminals, 
i.e. U([3, n/2)-U(0,O), U(f3,n)-U(p,n/2), U(f3, 3tt/2)- 
U(f3,ir) and U(f3, 0)-U(f3, 3tt/2), and then sum up the 
series in Eq. (|A2|) for a sufficient number of terms. 
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